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Abstract 

Threshold logarithms become dominant in partonic cross sections when the selected final 
state forces gluon radiation to be soft or collinear. Such radiation factorizes at the level of 
scattering amplitudes, and this leads to the resummation of threshold logarithms which ap¬ 
pear at leading power in the threshold variable. In this paper, we consider the extension of 
this factorization to include effects suppressed by a single power of the threshold variable. 
Building upon the Low-Burnett-Kroll-Del Duca (LBKD) theorem, we propose a decomposi¬ 
tion of radiative amplitudes into universal building blocks, which contain all effects ultimately 
responsible for next-to-leading-power (NLP) threshold logarithms in hadronic cross sections 
for electroweak annihilation processes. In particular, we provide a NLO evaluation of the 
radiative jet function, responsible for the interference of next-to-soft and collinear effects 
in these cross sections. As a test, using our expression for the amplitude, we reproduce all 
abelian-like NLP threshold logarithms in the NNLO Drell-Yan cross section, including the 
interplay of real and virtual emissions. Our results are a significant step towards developing 
a generally applicable resummation formalism for NLP threshold effects, and illustrate the 
breakdown of next-to-soft theorems for gauge theory amplitudes at loop level. 



1 Introduction 


It has long been known that the emission of soft and collinear gluons causes large corrections 
in perturbation theory, limiting its range of applicability. At the level of partonic cross sections, 
these corrections take the form of threshold logarithms of ratios of physical momentum scales, 
which diverge as the total energy or transverse momentum of emitted gluons (for soft or collinear 
emission respectively) becomes vanishingly small. Picking a dimensionless variable £ to measure 
the kinematic deviation from the threshold region, a differential cross section in this variable has 
a perturbative expansion of the form 
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We refer to the leading contributions in this expansion, characterized by the coefficients Cnm\ 
as leading-power (LP) threshold logarithms, while the Cnm specify next-to-leading-power (NLP) 
threshold logarithms, and the ellipsis denotes terms which are suppressed by additional powers 
of £. A variety of approaches for describing LP threshold effects exists in the literature, such 
as diagrammatic techniques based on factorization theorems m, approximations using Wilson 
lines 12 0 , renormalization group arguments [6], dedicated effective field theories mm and 
path integral techniques m Crucial to all these approaches is the notion of factorization, corre¬ 
sponding to the intuitive idea that long-distance effects in perturbative cross sections must have 
a universal nature, tied with the wave functions of the scattering states, and it must be possible 
to disentangle them from the contributions of short-distance high-energy exchanges. 

At the level of scattering amplitudes, long-distance effects at leading power manifest them¬ 
selves in the form of infrared and collinear divergences, which are known to factorize |121113) . For 
a generic n-point gauge theory scattering amplitude, this factorization takes the schematic form 
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Here the jet functions Ji (one for each external parton) contain universal collinear singularities 
which depend only on the color and spin quantum numbers of the external states, the soft 
function S ( n ) collects all soft singularities, which correlate all partons but do not depend on 
their energies and spins, and is a hard function, which is process dependent and infrared 
finite. One must finally divide each jet by its eikonal counterpart Ji, to correct for the fact that 
soft and collinear divergences have been double-counted. The factorization in Eq. (1.2) leads to 
exponentiation of soft and collinear factors in terms of a restricted set of anomalous dimensions: 
color singlet anomalous dimensions are known to three loops, since they can be extracted from 
the poles of partonic form factors O US]; the soft anomalous dimension matrix, on the other 
hand, is known at two loops for both massless and massive partons |16H25) ; the development 
of the necessary techniques to extend these results to the three-loop order and beyond is under 
way [26H39] . 

The factorization of leading-power infrared enhancements at amplitude level, embodied by 
Eq. (1.2), ultimately leads to the resummation of LP threshold logarithms for a large class 
of infrared-safe observables, under mild assumptions concerning the behavior of the associated 
real radiation m- In most cases, current knowledge of the anomalous dimensions allows one 
to perform this resunnnation up to N 2 LL or approximate N 3 LL accuracy. Given this increasing 
theoretical precision, and given the growing demands of current collider experiments, it is natural 
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to attempt to extend existing results to the next set of contribution^ in Eq. (1.1): indeed, these 
effects are known to be potentially significant [43W47] . as might be expected from the fact that 
they are still singular as £ —>• 0, albeit integrably so. Furthermore, there are theoretical reasons 
to expect that it should be possible to organize NLP threshold effects. Indeed, it has been known 
for a long time that soft gluon radiation effects at NLP can still be expressed in a universal 
way in terms of the non-radiative matrix element: this is the content of the Low-Burnett-Kroll 
theorem |48l 49] . This theorem was extended to the case of massless theories, where collinear 
effects become important, by Del Duca [50]. 

With these motivations, in recent years several steps were taken towards a deeper under¬ 
standing of NLP threshold effects and ultimately towards a resummation of NLP threshold 
logarithms [ID EMI, including results for a complete resummation of certain towers of NLP 
logs in specific processes. A fully general resummation prescription, applicable to all processes 
for which threshold logarithms can be resummed, has however yet to be developed. In particular, 
one must carefully disentangle the interplay between the soft expansion, which is well under 
control, and hard collinear effects, which can also be the source of threshold logarithms. 

In this paper, we will build on earlier results of Refs. HD ED ED- The first of these papers used 
path integral methods to perform a systematic expansion in powers of the momentum of emitted 
gluons, thus deriving a generalization of the soft-collinear factorization formula of Eq. (1.2), 
valid up to next-to-soft order. The second paper rederived these results using a diagrammatic 
approach, and also performed a first check of the formalism, applying it to the case of Drell-Yan 
production of a color-singlet vector boson, ft was shown that the next-to-soft Feynman rules 
correctly reproduce the known NLP logarithms in the double-real-emission contribution to the 
K factor for this process, up to NNLO [6al[69|. Whilst Refs. [HD 157] addressed the factorization 
and exponentiation properties of next-to-soft effects, the results were limited, as is the earlier 
work of Refs. [48L 49], by the fact that collinear singularities were not properly accounted for. 
In the case of Drell-Yan production, final state gluons are forced to be (next-to-) soft, so the 
interplay of soft and collinear enhancements arises only due to the presence of virtual hard- 
collinear gluons, which distort the emission spectrum of soft gluons at next-to-leading power in 
the soft expansion. In essence, in the context of Eq. (1.2), one must also include the possibility 
of soft emissions from the jet functions which contain all collinear radiation. 

The tools to tackle this problem, and thus to extend the Low-Burnett-Kroll theorem |481 4SJ 
to collinear-singular amplitudes, were developed by Del Duca [50] . who proposed a generalization 
of the standard hard-soft factorization, designed to organize the effects of next-to-soft radiation 
in the massless limit. A further step forward was recently taken by some of us in [64], where 
NLP threshold corrections to NNLO Drell-Yan production were precisely classified according to 
their soft or collinear origin, using the method of regions developed in Refs. [70lf72l J^] The results 
of Ref. [64] imply that in the case of electroweak annihilation cross sections the configurations 
needed to bridge the gap between the soft expansion and the threshold expansion are precisely 
the collinear-enhanced ones which are taken into account in Ref. [50] , 

Our aim in the present paper is to build on the results of Refs. HD ED ED EU, in the light 
of Eq. (1.2), to construct a precise formulation accounting for all NLP threshold logarithms for 
electroweak annihilation cross sections. With this in mind, we will improve upon the approach 


1 Recall that another important set of terms can be controlled to all orders, and can be shown to exponentiate 
in some cases: these are the contributions that have support only on the threshold, proportional to <5(£) if £ is the 
threshold variable mm- These terms will not be discussed here. 

2 The method of regions has also been recently applied, in a slightly different context, to characterise threshold 
effects in Higgs production via gluon fusion up to N'TO, in Refs. [731 174] . 


2 







of j50] by taking into account the factorization of soft modes (which were included in the hard 
interaction in Ref. m), by explicitly taking care of the double counting of soft-collinear sin¬ 
gularities, and by discussing the role played by the factorization vectors in the factorization of 
collinear from soft modes. This will lead us to a general expression for the radiative scattering 
amplitude, similar to that derived in [50], which contains all universal ingredients needed to or¬ 
ganize NLP threshold logarithms in the process at hand. We will then test our formalism, once 
again in the context of Drell-Yan production, by reproducing the known NLP logarithms in the 
mixed real-virtual contribution to the K factor for this process at NNLO, focussing for simplicity 
on the abelian-like term proportional to C\. We will describe how our formalism can be gener¬ 
alized to include purely non-abelian contributions, but we leave more detailed phenomenological 
analysis to future work. 

In parallel to the above developments, a recent body of work has examined the next-to-soft 
behaviour of gauge and gravity scattering amplitudes from a more formal point of view [75H95] . 
leading to a general theorem for the structure of next-to-soft corrections at tree-level and beyond. 
Particularly relevant in this regard is the study of Ref. [96], which approaches the classification 
of next-to-soft emissions from the viewpoint of soft-collinear effective theory (SCET), including 
also a detailed discussion of collinear effects. The relationship between this body of work and 
the approach of this paper has been recently clarified in Ref. m (see also [95]). Here we will 
follow up on this, by pointing out the implications of our results for loop corrections to next- 
to-soft theorems: our results broadly agree with, and substantiate with a concrete example, the 
arguments of Refs. [89, [96], showing that in the presence of collinear singularities next-to-soft 
theorems for scattering amplitudes require non-trivial corrections at one loop and beyond. 

The structure of our paper is as follows. In Section [2] we review in detail the modified fac¬ 
torization formula for next-to-soft emissions in the presence of collinear singularities proposed in 
Ref- [ED]) which we refine in several respects to emphasize the role of soft factors. As in Ref. [50], 
our final expression for the radiative amplitude at next-to-soft level contains a universal function 
describing the emission of a soft gluon from inside a jet, which we call the radiative jet func¬ 
tion. We calculate this function explicitly at one loop in Section [3] In Section [4] we assemble all 
necessary ingredients in order to reproduce the NLP threshold contributions to the Drell-Yan 
K- factor, providing a strong validation of our expression for the radiative amplitude. We discuss 
our results in Section [5] before concluding. Technical details are presented in an Appendix. 

2 Organizing NLP threshold logarithms 

In this section, we discuss in detail the construction of an expression, at the amplitude level, 
that directly leads to organizing threshold logarithms up to next-to-leading power. As already 
mentioned in the Introduction, a class of NLP threshold corrections arising from soft gluons 
at next-to-leading power in the soft expansion was analyzed in Refs. HU E23- The resulting 
factorization can describe those processes in which all threshold contributions arise only from 
soft emissions, but it will fail when collinear singularities are present in addition. This problem 
is shared by earlier works on next-to-soft corrections, such as the well-known investigation by 
Low [45]) who considered the emission of soft photons from hard scalar emitting particles, and 
showed how the resulting matrix element can be expressed in terms of (derivatives of) the non- 
radiative amplitud^] This seminal theorem was generalized to fermionic emitters in Ref. [ 49] , 
and the combined result is known as the Low-Burnett-Kroll theorem. This is still not the full 
3 See, for example, Appendix C of Ref. El for a modern review of this material. 


3 



story, however, as the results of Refs. |48l [49] apply only to particles with a non-vanishing mass 
m, so that one may take the energy E of soft gluons to zero in the parametrically well-defined 
limit E/m —> 0. One must clearly amend this argument for massless external particles, which 
in practice involves carefully disentangling collinear singularities (associated with m —> 0) from 
those associated with the soft expansion. This important generalization was achieved by Del 
Duca [50], who derived an expression for the radiative amplitude valid in the extended region 
m 2 /Q < E < m, where Q is the energy scale associated with the hard interaction. Taking 


m/Q parametrically to zero, one obtains a modification of Eq. (1.2), which is correct (at least in 


principle) up to next-to-soft level. For the purposes of studying threshold corrections, all these 
results are necessary, and in this paper we refer to this body of knowledge as the Low-Burnett- 
Kroll-Del Duca (LBKD) theorem^ 


Let us begin by reviewing the results of Ref. [50] . and placing them in the context of Eq. (1.2), 


which was developed in subsequent years. This will enable us to refine the original treatment, in 
particular explicitly extracting soft effects, studying the dependence on reference vectors defining 
the collinear region, and addressing the double counting of contributions which are both soft 
and collinear. We start our analysis by examining in more detail the soft-collinear factorization 


formula of Eq. (1.2) for the specific case of two hard colored particles, which we may write as 


A ( j^,a s (p 2 ),e 
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Ji (/3i,m,a s {fJ. 2 ),e) 
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Here we work in d = 4 — 2e dimensions, and pi (/?j) is the four-momentum (four-velocity) of 
the z-th hard particle. Again, Q is the energy scale associated with the hard interaction: for 
definiteness, we take here Q 2 = (p± + P2) 2 > 0. 


The partonic jet functions appearing in Eq. (2.1) are defined by m 


J(p,n,a s (n 2 ),e)u(p) = (0 |4> n (oo, 0)z/’(0)| p) , 


( 2 . 2 ) 


where i/j(x) is a quantum held inserted to absorb the external incoming parton with momentum 
p. The factor <f>„, on the other hand, represents a Wilson line stretching from the absorption 
point to infinity, along a direction fixed by an auxiliary vector nT, according to the definition 


< I , n(A 2 ,Ai) = V exp 


r A2 


ig s / dXn • A(Xn) 


' Ai 


(2.3) 


The presence of the Wilson line ensures that the definition is gauge-invariant. As was remarked 


in the Introduction, one must further introduce in Eq. (2.1) the eikonal counterpart of the jet 


functions, J ], in order to avoid the double counting of soft-collinear contributions. The eikonal 
jet is defined as 

J {ft, n, a s (p 2 ),ej = (0 |$ n (oo, 0)$g(0, - 00 ) | 0) . (2.4) 


Note that Eqs. (2.2) and (2.4) apply to incoming particles: for outgoing particles, one would have 


to reverse the direction of the Wilson lines. 

A key point of this discussion is the fact that one must associate with each leg of momentum 
Pi an auxiliary four-vector m. This ensures gauge invariance, but, more interestingly, it can 


4 The classification of NLP threshold effects from collinear regions has also been extensively discussed using 
effective field theory methods in the recent study [55] , 
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be physically interpreted as providing a means for measuring collinearity with respect to pi. 
The nfs are ‘factorization vectors’, and the full amplitude cannot depend on them, much as it 
cannot depend on the factorization and renormalization scale p. In fact, singular dependence on 
rii cancels between each jet and its eikonal counterpart, while non-singular n^-dependent terms 
cancel between the jets and the hard function. 

In order not to introduce in the jet functions Jj spurious collinear singularities not associated 
with emissions from the z-th hard parton, it is customary in factorization studies m to take 
rif 0. This has the advantage of allowing a cleaner identification of each subset of singular 
contributions, but it has the drawback of introducing a more complicated functional dependence 
on the factorization vectors. In contrast, a more physical viewpoint is to consider the rz, vectors 
as standing in to replace the other hard partons in the process: this is what is typically done in 
effective field theory calculations [96], and was also the approach followed in Ref. [64]. This choice 
leads to much simpler expressions, at the price of a degree of ambiguity in the classification of 
singular regions. Here, we will follow Ref. |50| . and define two dimensionless light-like vectors hj 
in directions opposite to pi, and such that 

nj = 0, fn-pi = Q. (2.5) 


In the calculations of Sections [3] and [4] it will in fact be more convenient to work with the 
dimensionful vectors 


. _ Q .. 

Ttn — 77 "i . 
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Hi -p 



( 2 . 6 ) 


note that for a two-parton process, such as that of Eq. (2.1), the definitions of Eq. (2.6) can be 
implemented simply by choosing 


n\ = p 2 , n 2 = pi, (2.7) 

an identification that will be useful in Section [4] 

Let us now turn to the soft function, which collects infrared singularities associated with the 
leading (eikonal) term in the momentum expansion of emitted gluons. For two-parton annihila¬ 
tion, it is defined by 


•S {fix • fi 2 ,a s {p 2 ),e) = (0|4> /32 (oo,0)4> /3i (0, — oo)|0) . (2.8) 


These expressions for the soft and jet functions implicitly define also the process-dependent hard 
function: for a given process, one defines T~L via Eq. (2.1), by taking all perturbative contributions 
to the soft and jet functions over to the left-hand side. The hard function then depends on the 
auxiliary vectors rij precisely in such a way as to cancel the finite parts of their contribution to 
the jet functions, so that the physical amplitude is independent of each m. 

As assembled in Eq. (2.1), the double counting of soft-collinear regions is solved by assigning 
soft-collinear poles to the soft function, and constructing ratios of jets and eikonal jets which 
contain only hard collinear singularities. An alternative arrangement of the same expression is 
to assign the soft-collinear poles to the jets, and construct the ratio of the soft function and the 
eikonal jets, which contains only soft wide-angle radiation. This ratio defines the reduced soft 
function, 

fix-fa _ i.. 2 \ \ S(fii- fi 2 ,a s {p 2 ),e) 


fix -nifi 2 - n 2 


, oi s (/i 


e = 


n iJiifii ■ ni,a s {p 2 ),e) ' 


(2.9) 


On the left-hand side, we have indicated that the arguments of the soft and jet functions must 
combine to form a ratio constructed to be invariant under the rescalings fii —> \fii- As discussed 
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Figure 1: Schematic depiction of the factorization of the amplitude into the non-collinear function 
H of Eq. (2.11) and external jet functions: (a) portrays the non-radiative amplitude, while (b) 
and (c) contribute to the radiation of an extra gluon. 


in detail in Haim El], for light-like /3j this invariance is broken for the soft function alone, as 
well as for the eikonal jets, due to the presence of collinear divergences in either factor. When 


the individual factors are combined into the reduced soft function, as in Eq. (2.9), collinear 


poles cancel and the invariance is restored. If, on the other hand, we work with light-like n*, the 
spurious collinear divergences associated with the Wilson lines in the rii directions do not cancel 
in S, so the expected invariance under the rescalings n* —> KiUi is not restored, as seen from the 


argument in Eq. (2.9). 


Making use of Eq. (2.9), we may now rewrite schematically the amplitude in Eq. (2.1) as 


A = T-L x S xJJjj, 


( 2 . 10 ) 


1=1 


where the functions {Jj} contain all relevant information associated with the collinear regions. 
Furthermore, in the remainder of this section, we will follow Ref. |50] and define a ‘non-collinear’ 
factor 


H = nxs, 


( 2 . 11 ) 


where the reduced soft function is absorbed into the hard function. The factorized structure of 
the amplitude is then as shown in Fig. [lja). Let us now describe how to generalise Eq. (2.1) to 
NLP level, building on Ref. m- First of all we wish to isolate the contributions to the radiative 
amplitude where the extra gluon is emitted by a collinearly enhanced configuration. With this in 
mind, and denoting the amplitude with an additional gluon emission by A one may naturally 
write 


A^k) = Ale^k) + A" e»{k) 


( 2 . 12 ) 


where we are suppressing color indices, is the polarization vector of the extra gluon, and Aj L 

(Ajl) represent emissions from the jet (hard) functions, respectively. The amplitude for emission 
from collinear configurations can be defined as 


2 2 
A'l = ^2 H{Vi — k\ Pj > nj) k, rii) Y[j{pj,nj) = ^Aft . (2.13) 

i =1 j^i i=l 


Here, for brevity, we have not displayed the dependence on the coupling and on e; we have 
introduced in H the notation of separating with a semi-colon the ‘active’ momentum (here 
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the momentum of the radiating leg) from the other vectors appearing in the function; most 
importantly, we have introduced the radiative jet function, 


J^{p,n,k,a s {n 2 ),e)u{p) = d d y e l(p k) ' y (0 | oo) if(y) j M (0) | p) , 


(2.14) 


representing the emission of a gluon from inside a jet function. It is useful to compare this 


definition with the non-radiative jet function defined in Eq. (2.2). Here the incoming particle (a 


quark of momentum p in this case) is absorbed by the held at position y, but, along the way, the 
insertion of the chromo-electric current j tl causes the loss of momentum k through the emission 
of a gluon. The gauge phase is then translated to infinite distance by the customary Wilson 
line. Finally the fact that fields are evaluated at different space-time points requires a Fourier 
transform in order to get the momentum-space correlator. For the purposes of the present paper 
we can take for the current the QED-like expression 


Ja(x) = ^{x)^T a if{x) 


(2.15) 


where T a is a color generator in the fundamental representation. In fact, we will continue to 
suppress color indices and matrices throughout the paper, except when especially relevant, since 
we will explicitly consider only QED-like diagrams and focus on contributions to the cross section 


proportional to Cp. Ultimately, however, we will need to generalize the definition in Eq. (2.14) 
to the full non-abelian theory, and to gluon-originated jets as well. A natural generalization of 


Eq. (2.15) is to consider the conserved (though not covariantly conserved) non-abelian current 
(see e.g. Ref. m) 

)a( x ) = fa bCF r( X ) A ^c(x)+'lf(x)^T a 'lf(x), (2.16) 

where f a bc are SU(N) structure constants. This choice has all the physical characteristics required 
for our definition: for example it allows radiation from a gluon jet, and it includes the proper 
non-abelian corrections to gluon radiation from a quark jet. There are, however, significant 
differences between the two cases concerning the application of relevant Ward identities and 
the renormalization properties, which modify to some extent the reasoning given below in the 
‘abelian’ case, without disrupting the general structure of the argument; therefore, we leave 
the full treatment of the non-abelian radiative jet to future work. We note that in either case 


Eq. (2.14) defines a universal object, depending solely on the properties of the given external 


particle (including its spin), but not on the details of the specific hard interaction process. We 
will explore further how to calculate this object in Section [3j here we begin by noting that the 


radiative jet function defined in Eq. (2.14), with the current in Eq. (2.15), obeys the simple Ward 
identity f50] 

Jf, (p,n, k,a s (p 2 ),e) = q J (p,n,a s (p 2 ),e) , (2.17) 

where q = ±1 according to whether the momentum p is incoming or outgoing respectively (q 
would represent the electric charge in a QED calculation). We may similarly consider the Ward 
identity for the entire amplitude: in this case, it takes the form 


WAp = 0 , 


(2.18) 


which in turn implies 

k y A* = -WA^. (2.19) 

We can use this result to relate emissions from inside the hard function, collected in the function 
Ajf, to those from the jet functions, represented by Aj r Indeed, note that so far Ajf has been 
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defined as a matching condition: the Ward identity shows that it is not independent from A^. 

( 2 . 20 ) 


To be explicit, we may use Eqs. (2.13) and (2.17) to find that 

2 2 


k^Afai, k) = ^2 Qi H (.Pi ~ k ',Pj, n j) II J (Pj, n j ) 
i= 1 j =1 

2 


2=1 


d 

H(pi',Pj,rij) + k/ 1 ( Qj^ H (Pi ~ k iPji n j) 


k-> o 


IPfo, 


n 


'31 i 


where we have Taylor expanded up to next-to-soft order in k in the second line. We may also 
rewrite the derivative in the second line as 


d 

dk^ 


H(pi - k-,pj,nj ) 


fc->0 


d 

dPi 


fi H (pi ; pj , rij). 


( 2 . 21 ) 


The zeroth order term in Eq. (2.20) now can be seen to vanish due to charge conservation (or 
rather color conservation in the QCD case), and the Ward identity of Eq. (2.19) implies 

•A%(pi,k) = '^q i (— TL H(p i ;p j ,n j ))Y[j(Pj,n j ). ( 2 . 22 ) 

i=i ^ op i ' j =1 


The reader may wonder at this point whether it might be possible to add a separate transverse 
contribution to the right-hand side. Such a contribution is argued to be absent at the level of the 
amplitude in Ref. [48] (see also [89, '90]), based on gauge invariance and locality considerations. 
Here we notice that the dependence of individual factors in Eq. (2.22) may in principle allow 
for transverse O(k 0 ) contributions, for example of the form k ,J ’/(rq • k). All such contributions 
must however cancel in the complete amplitude, and we will not include them in our analysis. 

Returning now to the amplitude for emission from the jet functions, it is convenient, following 
Refs. |50( dOlj . to introduce, for each jet, a decomposition of the polarization smrj^jas 


rf v = + K 111 ', 




(2.23) 


so that G^ u satisfies 


P IJ G IW = 0(k) , G/jyk 1 ' = 0. (2.24) 

This decomposition leads to a particularly simple expression for the emission of a ii-gluon from 
the jets. Indeed, the Ward identity in Eq. (2.17), together with Eq. (2.13), yields 

K T = G k \,2 H(Pi~k] pj ,nj)f[j ( pj , nj) 


= Qi 


2pi ■ k — k 2 

(2 pj - kY 
2 pi - k — k 2 


o =i 


d 


A - ( nfa’Pj* n j)) n j (p 


uj , rij 


3=1 


(2.25) 


where K 1 - 1 ' 1 is the K tensor appropriate to the i-th jet, we have again Taylor expanded in k, and 
we have recognized the non-radiative amplitude A in the second line. Upon combining this result 


with the emission from the hard function, as given in Eq. (2.22), one obtains 
2 2 


A H ^ + Y J ^'K7' = 

2=1 


2=1 


(2 Pi - kr 

2 pi - k — k 2 


" 4+ ( G 7-?L H (PiiP3i n 3) ) Ip(W.»i) 


dp[ 


3 


. (2.26) 


5 Our definition of the G and K tensors differs from that of Ref. !50] due to the fact that we have taken the 
momentum k to be outgoing. Note also that in our calculation k 2 = 0, although most of the argument goes 
through also for off-shell k 
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Using Eq. (1.2), one may rewrite this in a more useful form, as 


A H ^ + ^ A'l 


KV L 


2=1 




( ( 2 P* ” k Y | 9 \ a 

\2 Pi -k-tf~ ' dp" j 

d 2 

- H{p i -p j ,n j ) G U A — IJ J(Pj,nj) , 


(2.27) 


so that derivatives with respect to hard momenta act only on the full non-radiative amplitude, 
or on the process-independent jet functions. 

We must now, finally, consider the emission of a G-gluon from the jets, which is given by 


2 2 

*22 A J u l G V A = ^2 G T H C Pi ~ *5 Pj’ n j) J APii k , n i ) J! J ( Pj, 


n 


'0 ) ■ 


(2.28) 


2= 1 


2=1 


As argued in Ref. |50] . projection with the G tensor ensures that this contribution starts at 
next-to-leading power in the soft expansion. One may therefore Taylor expand the shifted hard 
function, keeping only the zeroth order term. Adding this to Eq. (2.27), the radiative amplitude 
finally becomes 


A^(pj,k) = J2 


2=1 


a G ( 


(2.29) 


+ n(pj , n 3 ) S(/Jj , n 3 ) CfA | ® -J^pJ(pi,ni) ) JJ J(jPj, nj) 

Pi 7 


where we have restored the factorized expression for the non-collinear function H in terms of the 
hard factor % and the reduced soft factor S. 


In the remainder of the paper, we will provide strong evidence that Eq. (2.29), first derived 


with slight modifications in Ref. [50], provides the missing ingredient for a complete reconstruc¬ 
tion of NLP threshold logarithms in electroweak annihilation cross sections, generalizing the 


classic leading-power factorization formula given in Eq. (1.2). Some comments are in order. 


We expect, and we will confirm below, that when Eq. (2.29) is used to construct a NNLO 


annihilation cross section, by contracting it with the corresponding tree-level matrix ele¬ 
ment and integrating over phase space, it will correctly generate all threshold logarithms 
up to NLP. Phase space integration does not present difficulties, since next-to-soft contri¬ 


butions to Eq. (2.29) can be integrated with the well-known factorized expression for the 


phase-space measure at leading power, while correction to phase space are only needed for 
leading-power matrix elements, and were discussed in detail in Ref. m 


• A new universal quantity arises at this level, the radiative jet function, defined originally 
in Ref. |50j. In order to predict NLP logarithms in specific scattering processes, one must 
calculate this quantity, which is the subject of the following section. 


The only process-dependent contributions on the right-hand side of Eq. (2.29) are the full 
non-radiative amplitude A and the hard function PL. 


• At the level of singular contributions, our result complements the analysis of Ref. [50] by 
correctly accounting for the soft and soft-collinear regimes, and in particular by subtracting 
the double counting of the latter with the inclusion of eikonal jet functions. 
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It is important to note that our treatment suffices for electroweak annihilation cross sec¬ 
tions, which do not involve final-state QCD jets, since all radiated gluons must be (next-to-) 
soft in this case. In the presence of final-state hadrons (for example in the case of Deep 
Inelastic Scattering (DIS), NLP threshold logarithms may in principle be associated with 


hard collinear emission. These are potentially not taken into account in Eq. (2.29), which 
relies upon the soft expansion. 

Having now presented our general framework, we turn to the calculation of the radiative jet 


function defined by Eq. (2.14), focusing in the present case on the abelian-like contributions 


generated by the current in Eq. (2.15), up to one loop. 


3 The radiative jet function 


Eq. (2.29) expresses the radiative amplitude for electroweak annihilation in terms of a process- 


dependent hard function, together with a number of universal quantities. The soft function, the 
non-radiative jet and its eikonal counterpart are already defined at leading power in the soft 
expansion, and can in principle be taken over without changes. One also encounters, however, 


a new universal spin-dependent quantity: the radiative jet function of Eq. (2.14). This quantity 


is defined for finite emitted energy k, and must be computed with NLP accuracy in the soft 
expansion, while keeping control of all collinear singularities which might interfere with the soft 
emission: all these contributions are expected to give rise to threshold logarithms when contracted 
with the tree-level complex-conjugate amplitude and integrated over the soft gluon phase space. 

To substantiate our claims, in this section we will compute the radiative jet function to one 
loop order, selecting all terms that will contribute abelian-like threshold logarithms, proportional 


to Cp, to the NNLO cross-section. We will therefore choose the current in Eq. (2.15), omit graphs 


containing three-gluon vertices, and pick the appropriate color structure for vertex corrections 
which would contribute terms proportional to CaCf, which we disregard. Before turning to that 
calculation, we will however pause to briefly discuss and motivate our choice of reference vectors 
nT, which are necessary ingredients for non-radiative, eikonal and radiative jets. 


3.1 Jet functions for light-like 


It is customary in the factorization literature to compute the jet functions appearing in Eq. (1.2) 
with a reference vector off the light cone, n 2 ^ 0 (see, for example, Ref. H2I). This has several ad¬ 
vantages, in particular in the case of all-order proofs and for renormalization group studies. Most 
notably, keeping nT off the light cone prevents the appearance of collinear divergences associated 
with the nT Wilson fine, which are spurious and should in general be removed; furthermore, 
keeping n generic allows to test perturbative calculations by verifying their independence on n; 
finally, space-like n preserves certain analytic properties of the amplitude which are useful for all¬ 
order analyses. The non-radiative jet function J, and its eikonal counterpart J % are well known 
and they have been computed at one loop for arbitrary (non-null) n for example in Ref [I, 13). 

Here we would like to argue that in the present case there are considerable computational 
advantages to be gained by keeping n 2 = 0, and at the same time the issues of interpretation, that 
can be quite subtle for generic correlators, can easily be dealt with for the radiative amplitude 
that we are discussing here. To illustrate the problem, let us begin by considering the one-loop 
calculation of the non-radiative jet function J(p,n). For an external incoming quark, one must 
consider the diagrams in Fig. |2j where the second graph denotes a UV counterterm associated 
with the renormalization of the vertex, and we have omitted the similar counterterm graphs 
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Figure 2: Feynman diagrams contributing to the one-loop jet function. Here Jy^ CT denotes the 
counterterm associated with the vertex graph, Jy\ while counterterms associated with external 
leg corrections have been omitted. 


associated with external leg corrections. For n 2 7 ^ 0, the complete jet function must depend 
on the dimensionless variable z = (jp ■ n) 2 /(n 2 p 2 ), as a consequence of the invariance of the 
eikonal Feynman rules under the rescaling n^ —> nn^. Furthermore, the function is well defined 
in dimensional regularization due to the presence of the energy scale p ■ n. For n 2 = 0, both of 
these properties are lost: dependence on p ■ n is in principle ruled out by rescaling invariance, so 
that all integrals arising in the relevant Feynman diagrams are effectively scale-less, and must be 
defined to vanish in dimensional regularization. In fact, there is an extra twist: while the diagrams 
vanish, one finds that there is a residual dependence on p ■ n in the UV counterterms, due to an 
anomalous breaking of rescaling invariance originating from the collinear pole associated with 
emission from the light-like Wilson line [151 IT9j . 

To give a concrete example, consider the first diagram in Fig. [ 2 J which we denote by Jy\ In 
d = 4 — 2e dimensions, and for n 2 = 0, it is given by 


Jy\p,ri;e) = 


2i P 2e g 2 s 


2i/i 2t g 2 s 


d d k (jl- ff) ji 

(2ir) d k 2 2 n ■ k(p — k) 2 

d d k f 1 , f 1 , 
dx / dy 


(3.1) 


2 y (tf- ft) ft 


(2n) d J 0 J o [yfc 2 — 2xyk ■ p + 2(1 — y)n ■ k\' 


In the second line, we have introduced Feynman parameters, and we note the characteristic 
parameter dependence of the k 2 term in the denominator, which arises in the presence of linear 
denominators. Carrying out the momentum integration and using the Dirac equation, one may 
rewrite this as 


Jy\p,n-,e) = ^ ( 47 r/i 2 ) e T(1 + e) (~2p 


n) 


1 


dyy- l+e {l-y) 


—l—e 


(3.2) 


At this point, one might be tempted to interpret directly the y integral as B(e, — e) = 0. More 
accurately, one observes that the integral is not well defined for any values of e, and must 
therefore be defined to vanish in dimensional regularization. In this simple case, it is actually 
easy to disentangle the ultraviolet divergence (arising from the region y —> 0 ) from the infrared 
one (arising from the region y —> 1). One may simply insert a factor of (1 — y) + y = 1 to 
see explicitly that infrared and ultraviolet poles cancel exactly; using different regulators yields 
compatible results. The standard treatment at this point is compute the renormalized jet function, 
isolating the ultraviolet divergence in Eq. (3.2) and defining the appropriate MS counterterm to 


subtract it. The sum of Jy ^ and Jy’ CT is then simply the negative of the UV pole of Jy‘ 
which is correctly interpreted as an infrared divergence. For n 2 = 0 one recovers the expected 
double soft-collinear pole because of the explicit collinear divergence in Eq. (3.2). The external 
leg corrections in Fig. [2] can be treated similarly. In what follows, we will take the alternative 


A i) 


r( l) 
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route of computing bare jet functions, which will allow us to simplify considerably all calculations, 
and make a more direct contact with the calculation performed with the method of regions in 
Ref. |64] , and eventually with the SCET approach described in Ref. [96]. This approach works for 
our current problem, and for the choice n 2 = 0, because of the simple renormalization properties 


of Eq. (2.29). More specifically, we note that the left-hand side of Eq. (2.29) is an on-shell 


scattering amplitude, and thus it is not renormalized: all counterterms needed for the various 
factors must cancel. Furthermore, the right-hand side can be expressed almost entirely in terms 
of the complete non-radiative amplitude A, which is also renormalization-group (RG) invariant. 


The first line of Eq. (2.29) is already in this form; the second line can be brought to a similar form 


simply by noting that, in each term in the sum over external partons, only one non-radiative jet 
function is missing to reconstruct the full non-radiative amplitude. Multiplying and dividing by 


that jet one can write Eq. (2.29) as 


k) = Y \ ® ' V ,o + GT 


i =1 


( ( 2 Pi - kY 
\ 2pi ■ k — k 2 


d 

dp\ 


(3.3) 


+ G7 


J u (pi,k,m ) d . 

~T( - 1 - G 1 nJ{pi,ni 

JyPii n i) dp7 


A(pf,Pj) ■ 


In this form, it is evident that the factor in square brackets in the second line must be RG 
invariant by itself, to all orders in perturbation theory. This is easily verified by inspection, 
noting from the diagrammatic expansion that the UV divergences of the radiative jet function 
Jfj, are the same as those of the non-radiative jet J, simply multiplied by the tree-level expression 


for Jn- We are thus free to compute Eq. (3.3) in either bare or renormalized perturbation theory. 
The advantage of using a light-like reference vector is now apparent: with m? = 0, radiative 
corrections to the bare non-radiative jet function vanish to all orders in perturbation theory, and 


one can simplify Eq. (3.3) by setting 


J{pi,rii ) = 1 . 


(3.4) 


The term containing the derivative of the non-radiative jet function now vanishes, and we are 
left with the simple expression |^] 


■A ll (Pj, k) = Y 


i =1 


Ql 2pf- k -l 2 + Qi + ) A te'P3) ■ 


(3.5) 


We will use this explicit form in Section [4] and we believe that the simple formal properties 


embodied in Eq. (3.3) and Eq. (3.5) will be useful also for the all-order analysis that will have 


to be performed in order to construct a resummation procedure for NLP threshold logarithms. 


3.2 The radiative jet function at one loop 


The main result of Section [3.1| is that, for the purposes of the present calculation, we are allowed 
to work with bare quantities, and with light-like reference vectors for jets. This has considerably 
simplified our task, since perturbative corrections to non-radiative (eikonal) jet functions can 
be taken to vanish to all orders. The same is of course not true for the radiative jet function 


Strictly speaking, one could write Eq. (3.51 explicitly in terms of the ratio = J ll (pi,k)/J(pi). As 


Eq. (3.41 makes clear, J ll (pi,k) = in dimensional regularisation. As well as being UV finite, the ratio is 

also free of soft (but not collinear) singularities due to the correspondence between UV and IR poles. If a different 
regularisation scheme is used, these properties persist for J M , but not for itself. 
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Figure 3: Contributions to the QED-like terms in the one-loop bare radiative jet function. 


defined in Eq. (2.14), which in fact depends on several momentum scales. We will now evaluate 


this function at one loop, using a light-like reference vector, and working in bare perturbation 


theory as required by our reasoning in Section 3.1 


We begin by defining the perturbative coefficients of the radiative jet function via[^] 

oo 

J u 0, n, k ; a s , e) = g s ^ ) J { u n) (p, n, k ; e) . (3.6) 


n =0 


The operator definition in Eq. (2.14) gives a straightforward result at tree level, consisting of a 


single emission from the external leg of momentum p, and yielding the expression 


2p ■ k p ■ k 


p 


+ 


k h 


i k n 


p ■ k 2 p ■ k 2 p ■ k 


(3.7) 


In the second line, we have chosen to decompose the result into spin-independent and spin- 
dependent parts, where the latter takes the form of a magnetic moment coupling to the fermion 
leg, involving the Lorentz generators 


= 


4 L 


7 > 7 


(3.8) 


The expression in the second line of Eq. (3.7) naturally generalizes to the case of emission from 


hard partons of different spin, simply by choosing the appropriate form for the corresponding 
Lorentz generator. 

At one loop, things are more complicated, and the relevant contributions are depicted in 
Fig. i Diagrams (a) and (b) have a more intricate kinematic dependence, as they link the physi¬ 
cal external leg to the Wilson line, while diagrams (c) and (d) are standard radiative corrections 
to the emission amplitude. External leg corrections are not included since they are pure coun¬ 
terterms. The diagrams linking the parton line and the Wilson line are readily evaluated. For 
our purposes, we need to include terms up to O(e), since phase-space integration over k will gen¬ 
erate singularities. On the other hand, we can neglect terms quadratic in k^ since they cannot 
contribute to NLP threshold logarithms. The result is 


M 1 ) 

J a+b 


(P,n,k;e) 


(2 p-r 


- + 4 + 8e 


n ■ k p u 
p ■ k p ■ n 


ftY _ n v 
2 p ■ k p ■ n 


7 Given that we are explicitly talking about QCD corrections from now on, we have replaced the electromagnetic 
charge factor qi in section [2] with the strong coupling g a . Note that in our convention the sign of q; is negative for 
an incoming quark line. 
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+ <i + 3 £ )(-^ + Yi_YlA 

p • k p • n p • fcp • n 


(3.9) 


Notice that this combination has only a single pole, while the individual diagrams (a) and (b) 
have (canceling) double poles. Diagrams (c) and (d) do not depend on the reference vector n, 
and they yield 


•4+d (p, n, k ; e) = (2 p ■ k)~ 


ftY k v 


p ■ k p ■ k 
+ e 


5 ftY k 

+ x ——— + 


2 p ■ k p ■ k 
k v 


(3.10) 


k'f 

5 ^T + 2 t 

p ■ k p ■ k 


The full result for the bare one-loop radiative jet function is then given by 


(p, n, k ; e) = (2 \p ■ k)~ 


- + 4 + 8e 


n ■ k p v 


n 


+ 


1 1 

-3e 

e 2 


p ■ k 


+ (1 + 3e) 


p ■ k p ■ n p ■ n 

Y p v ft i 1 


- (1 + 2e) 


\k a Y. av 
p ■ k 


p■n p■kp■n 


+ 0(e 2 , k) 


(3.11) 


where we have again chosen to write this in terms of spin-independent and spin-dependent 
components, for the n-independent terms. A comment is in order regarding the e pole appearing 
in this result. As explained in the previous section, we have avoided the introduction of UV 
counterterms given that these will cancel when all ingredients in eq. (3.5) are combined. It 
follows that the singularity appearing here can be interpreted directly as the infrared pole that 
will appear in the final amplitude. 

An important check of the above results is the Ward identity, Eq. (2.17). We note first that 
it is verified by the tree-level jet emission function in Eq. (|3.7|): indeed 

(3.12) 


k v J v(si) (jp,n,k) = -g s j(°\p,n ), 


where we have used the fact that, at tree-level, J^°\p,n) = 1, and that our expression has been 
derived for k 2 = 0, as well as our convention for the sign of the charge of an incoming fermion. 
At the one-loop level, the results of Eqs. (3.11) and 3.12) imply that 


k u J u ^ = 0 = - g s J {l) 


(3.13) 


given that loop corrections to the non-radiative jet function vanish. Note again that we are 
verifying the Ward identity only up to corrections ~ 0{k 2 ), since our jet functions have been 
computed for on-shell k. 

In Eq. (3.3), one must finally contract the jet emission function with the G tensor defined in 
Eq. (2.23). The reader may verify that, remarkably, the one-loop expression in Eq. (3.11) is an 
eigenstate of G. More precisely one finds 

G^jW (p } n , k) = (p, n, k) . (3.14) 

We now have everything we need in order to test our factorized expression for the radiative 
amplitude, Eq. (3.5), in the context of a NNLO calculation. To this end, we consider Drell-Yan 
production in the following section. 
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4 Application to Drell-Yan production 


In the previous sections, we have assembled the ingredients needed to organize threshold cor¬ 
rections at next-to-leading power, at the level of scattering amplitudes. As already mentioned 
above, we will now test our results in the concrete example of Drell-Yan production at NNLO, 
with the goal of reproducing abelian-like logarithms (weighted by the color factor C \) in graphs 
involving the real radiation of a (next-to-) soft gluon, dressed by a virtual correction. This is a 
non-trivial check that our formalism is correct, given that one expects an interplay between soft 
and collinear contributions when virtual gluons are included. 

Drell-Yan production is the simplest testing ground for our formalism, since it involves only 
two colored particles at leading order, it has a non-trivial abelian-like limit, and it does not involve 
final state jets; furthermore, NNLO corrections have been known for many years |45l [65] I102L 
1114]. A more interesting process for phenomenological applications is Higgs production via gluon 
fusion, where NLP logarithms are now known at N 3 LO [73, 74] (see also refs. [ 105 , 106 ]): we plan 
to tackle that process in future work, when the full non-abelian generalization of our formalism 
has been worked out. 

We recall that, for the purely real emission contributions, involving two real gluons at NNLO, 
all threshold logarithms at both LP and NLP have been shown to arise from the soft expansion, 
with no contamination from collinear singularities, as has been verified up to NNLO in Ref. m- 
Mixed real-virtual corrections at NNLO, which we examine here, are the first instance in which 
the collinear generalization of Low’s theorem is required. In a diagrammatic framework, the 
various contributions involved have been analysed recently in Ref. [64], using the method of 
regions of Refs. |70!f72| . We will make contact with this analysis in what follows. 


4.1 The real-virtual Drell-Yan K-factor at NNLO 


We begin by very briefly setting up or notations and conventions. The leading-order Drell-Yan 
process for the production of an off-shell vector boson V* of invariant mass Q 2 proceeds through 
the process 

q(p) + q(P) -> v*(Q), ( 4 . 1 ) 


where q and q denote a quark and antiquark respectively, and arguments label their four- 
momenta. Defining the squared centre-of-mass energy s = (p + p) 2 , one may introduce the 
dimensionless threshold variable 


z = 


91 

s 


(4.2) 


representing the fraction of available energy carried by the final state photon; the threshold limit 
then corresponds to z —> 1. The K-factor at fixed order in perturbation theory is defined by 


K (n \z) 


1 da^ n \z) 
er(°) dz 


(4.3) 


where is the n-loop Drell-Yan cross section. For the case of emission of a single real gluon 
with momentum k, one may introduce the Mandelstam invariants 


t = — 2k ■ p , u = — 2k ■ p , 
which may be parametrized as [68] 

t = - 2s(l - y)(l - z ), u = -2sy(l-z), 


(4.4) 


(4.5) 
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where 0 < y < 1. Using these variables, and setting, for simplicity, the renormalization scale 
fi 2 = Q 2 , the real-virtual contribution to the NNLO K-factor can be written as 


Kg\z) = 




dy[y(l-y)\ c JQ w A t + A\A 


(4.6) 


where A IV and A r are the amplitudes for single real gluon emission, at NLO and LO respectively. 

Although the complete NNLO Drell-Yan K-factor has been known for a long time [68] . no 
separate result exists in the literature for the real-virtual contribution. We have reproduced the 
relevant calculation, and we give here the resulij^] for threshold contributions up to NLP. One 
finds 


A'(?( 2 ) = (^C F 


32 


V 0 (z) - 1 


16 
H—v 


— 4X>i(z) + 3V 0 (z) + 4 L(z) - 6 


4 r 


16 V 2 (z) - 24Pi(z) + 32 V 0 (z) - 16 L 2 (z) + 52 L(z) - 49 


+ - 
e 

128 

- —V 3 (z) + 96 V 2 (z) - 256Pi (z) + 256P 0 (^) 

O 

1 OQ I 

+ —L 3 (z) - 232 L 2 (z) + 412 L(z) - 408 L 


(4.7) 


where we defined 


V n (z) = 


log n (l - z) 


L(z) = log(l - z). 


(4.8) 


(1 - z) 

Further details of the calculation can be found in Appendix [A] Note that we do not include 


the Dirac 5-function contributions in Eq. (4.7), as these mix with virtual corrections that are 
irrelevant for the purposes of this paper. 


4.2 Reconstructing the K-factor at NLP 


We are now going to reproduce Eq. (4.7), starting from our factorized expression for the ampli¬ 


tude, given in Eq. (3.5), and including the required contributions from the phase space measure 


at NLP in the threshold expansion. Since the calculation is not technically difficult, we sim¬ 
ply outline the main steps and give the results for the various contributions to the K-factor, 
distinguishing their physical origin. 


4.2.1 NLP corrections to the phase-space measure 


Before discussing how the T> n (z) and L n (z ) terms in (4.7) are reproduced in our factorized 
approach, we address the issue of how the amplitude factorization presented in Eq. (2.29) and 
Eq. (3.5) can be employed to obtain cross section result^ Eq. (4.6) expresses the required 
contribution to the K-factor in terms of the density 


V = A\ v A r + A* A IV 

= Plp + Pnlp + • • ■ , (4.9) 

8 Here and throughout, for brevity, we neglect terms involving transcendental constants, as was done in Refs. ( 57 1 
164) . Such terms do not contain any new information. 

9 Note that the present situation differs from that of Ref. [73] , where the three-loop results for the form factors 
contributing to Higgs and Drell-Yan production could be used to infer the coefficients of the distributions V„(z) 
at N 3 LO for the corresponding cross sections; in fact, that procedure does not extend to the L n (z) terms. 
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where the second line denotes the expansion in powers of 1 — 2 up to NLP. This must be combined 
with the overall phase-space prefactor 

z e (l-z) l ~ 2t = (1 -z) l ~ 2t 

which arises from the real gluon phase space. 

Working at next-to-leading power in (1 — z) has the advantage that corrections to the phase 
space measure and to the matrix element do not interfere. One may, in fact, write the differential 
cross section schematically as 


l + e(l-z) + 0({l-z) 2 ) 


(4.10) 


da = d<& 3i LP (T’lp + T’nlp) + ^^nlp^lp > 


(4.11) 


where d $3 (N)LP denotes the three-particle phase space at the required order in the threshold 
expansion. In the present context, we see that we can proceed by integrating the full squared 


matrix element, given by Eq. (3.5) contracted with the tree-level amplitude, with the eikonal 


expression for the phase-space measure; we can then include phase-space corrections, where 
however only the leading-power squared matrix element is needed. Note that NLP corrections 
to the phase space measure affect the cross section to all orders, starting with the tree-level 
emission of a single gluon, as was already noted in Ref. m in order to reproduce the known 
NLP logarithms in the one-loop Drell-Yan K-factor. Appendix B of that paper also included a 
discussion of how to obtain the NE contribution to a general nmlti-gluon phase space. In what 
follows, given the simplicity of the one-gluon-emission calculation, we will not make explicit use 


of Eq. (4.11), which however would be useful in order to construct all-order expressions. 


4.2.2 From tree level to one loop 

In order to reproduce all NLP terms in the real-virtual contribution to the NNLO Drell-Yan 


K-factor, we must calculate all contributions appearing in Eq. (3.3), before integrating over the 


real-gluon phase space according to Eq. (4.6). Before performing this calculation at one-loop 


level, it is worth pausing to briefly remark upon what happens at tree level. In that case, the 
derivatives in Eq. (3.3) and (3.5) act on momentum-independent leading-order quantiticsj^j and 
thus vanish. The factor including the radiative jet function is easily computed using Eq. (|2.23 ) 


and Eq. (3.7), and it gives 


G^4 0) (p,n,fc) = J (0) (p,n), 

P ■ 


(4.12) 


where we have used the fact that = 1. Using the leading-power factorization formula 


in Eq. (1.2), one then writes 


_ i k 

'H(pi,n i )S(pi,n i )G” fl J I ,(pi,ni,k)Y[j(pj,n j ) = - A(p t ,pj ). (4.13) 

Pi ' ft 

37 =* 


This combines with the first term on the right-hand side of Eq. (3.5) to reproduce the non 


radiative amplitude dressed by a complete (spin-dependent) gluon emission. Upon integrating 
over the real gluon phase space, this precisely reproduces the NLO K-factor calculation, already 
carried out using an effective approach in Ref. I5ZI- 


10 Note that the derivatives do not act on the momentum dependence implicit in the particle wave functions, 
for example the quark spinors in this case. 
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Returning now to the one-loop example, it is useful to distinguish three contributions to the 
cross section which are physically distinct. Let us begin by rewriting the one-loop contribution 


to Eq. (3.5) as 




2—1 




+ (Pi , k)A {0) ( Pi]Pj) 


(4.14) 


where we expanded the amplitudes in powers of a s /( 47 r), as in Eq. (3.6). Using Eq. (4.12), one 
finds that the first and third terms on the right-hand side combine to give the non-radiative 
amplitude, dressed by a complete (spin-dependent) gluon emission. This is directly analogous 
to the tree-level behaviour described above. Together with the remaining terms, there are then 
three contributions to calculate: the dressed non-radiative amplitude, the derivative of the non- 
radiative amplitude, and the radiative jet contribution. We compute them in turn. 


4.2.3 The dressed non-radiative amplitude 

The one-loop non-radiative amplitude can be easily obtained from the well-known result for the 
one-loop quark form factor (see, for example, Ref. [L3]). Setting p? = Q 2 and reabsorbing the 
factors of ln( 4 - 7 r) and 7 e in the renormalization scale, as appropriate to the MS scheme, it is 
given by 

= ~ £ r ( r(l- r 2e) +e) (| + ~ e + 8 + 166 + ° (e2) ) ’ (415) 

where the factor z e arises since the non-radiative amplitude must be computed with (p + p) 2 = s. 
As is implicit in the above discussion, this dresses the tree-level amplitude involving one real 
gluon emission. The relevant amplitudes were reviewed recently in Ref. m- The results for the 
squared matrix elements at leading and next-to-leading power are 

(4.16) 


|A,lp | 2 = 16(1 -e)g 2 s S ~, 


and 


•A",NLP “^r,LP + -^LP “4r,NLP — 8(1 e) g s ^ , 


(4.17) 


respectively, where again LP and NLP refer to the order in the expansion in powers of (1 — z). To 
obtain the appropriate contribution to the Drell-Yan K factor at NNLO, one must now multiply 


the sum of Eqs. (4.16 4.17) by the form factor contribution in Eq. (4.15). After carrying out the 


phase space integration for the emitted gluon according to Eq. (|4.6|), the result is 

V 0 (z) - 1 


tfdl’M = 


r^) 2 n r 


+ ~2 
e z l 


— %V\(z) + 6Vq(z) + 8 L{z) — 14 


16 


4 V 2 {z) - 6 P 1 (z) + 8 V 0 {z) - 4 L 2 (z) + UL(z) - 14 


+ — 

e 

128 

- —V 3 (z) + 96 V 2 {z) - 256Pi (z) + 256 V 0 (z) 

O 

19s I 

+ -yL 3 (z) - 224 L 2 (z) + 448 L(z) - 512 \ . 


(4.18) 


Not surprisingly, all plus distributions are correctly reproduced by the dressed Born amplitude: 
the derivative and radiative jet contributions are expected to be strictly NLP, as we verify below. 
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4.2.4 The derivative of the non-radiative amplitude 


The second term on the right-hand side of Eq. (3.5) consists of the non-radiative amplitude, 


differentiated in turn with respect to each external momentum. Considering for example the 
derivative with respect to p, Eq. (|4.15 ) gives 


G^(p,k) 


dp v 


—= -r + 
p-p 


p-k . 


p ■ k 


P 


A™ 


(4.19) 


This derivative, and the corresponding one with respect to p, must be contracted with the tree- 


level real-emission amplitude, before being integrated over phase space as in Eq. (4.6). Including 
also complex conjugate graphs, one finds that the contribution to the K factor is 


k qI{z) = l^C F ) + — — 4L(z) + 3 + 64:L z (z) — 96L(z) + 128 


O'.s 


4-7T 


32 16 


e l 


(4.20) 


which contributes at NLP, as expected. 


4.2.5 The radiative jet contribution 


The last part of the NLP K-factor comes from the second line of Eq. (2.29), and is due to 


those radiative jet contributions that have not yet been included in the dressed non-radiative 


amplitude. These can be extracted from Eq. (3.11). At this stage, it becomes inevitable to make 


a precise choice for the factorization vectors nf in the radiative jet functions. In the leading- 


power factorization formula, Eq. (1.2), it is possible to engineer a complete cancellation of the 


n dependence of the non-radiative jets: n-dependent poles cancel between jets and eikonal jets, 
while n-dependent non-singular terms are cancelled by the hard function PL. At NLP, we have 
already made a special choice by setting n? = 0; with this choice the presence of n-dependent 


poles in Eq. (3.11) suggests that in order to achieve formal independence from n we would need to 


introduce a further subtraction, replacing collinear poles in the direction m with those associated 
with the physical momentum of the parton colliding with parton i. There is however a simple 


shortcut, already introduced in Eq. (2.7): we may simply choose 


n = p, 


n = p. 


(4.21) 


This is physically motivated by the fact that the Wilson line in the direction m acts as a re¬ 
placement for the parton colliding with parton i; furthermore, considering for example emission 
from the parton with momentum p, one observes that p is the only other light-like vector in the 
problem; finally, the method-of-regions calculation of Ref. (64j shows that this choice of refer¬ 
ence vectors does indeed account for all the poles in the radiative amplitude. For processes with 
more legs, one would then choose rij to be the anti-collinear vector appearing in the Sudakov 


decomposition related to leg i , i.e. satisfying Eq. (2.5), with Q the appropriate hard scale. 


With the choice in Eq. (4.21), we proceed by contracting Eq. ( 3. 11|) (for both the p and p 
legs) with the NLO tree-level amplitude, including complex conjugate diagrams, and then we 
integrate over phase space. We find 


K 


( 2 ) 
collinear 


(*) = 



12 L(z) - 5 - 72 L 2 (z) + 60 L(z) - 24 


(4.22) 


where the label ‘collinear’ for this contribution will be discussed in the following section. Once 
again, as expected, this contribution is strictly NLP. 
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4.3 Discussion 


We have now calculated all necessary ingredients to reproduce the real-virtual, abelian-like con¬ 
tribution to the NNLO Drell-Yan K factor. Adding these together according to 


Kg’tz) 


Km(z) + 4%) + K‘H rme Jz), 


(4.23) 


and using the results for each term given in Eqs. (|4.18 4.20 4.22), we precisely reproduce the 


full perturbative result, given in Eq. (4.7). This is a non-trivial check of the validity of the NLP 


factorization formula, Eq. (3.5), particularly given the fact that all terms in the cross section 


(including even the constant in the finite part) are reproduced: no complete prediction for single- 
logarithmic contributions at NLP and NNLO has been given before. 

At this point, it is interesting to compare our results with the recent study of Ref. [64], which 
used the method of regions [MS] to classify the same contributions to the NNLO Drell-Yan 
K-factor. The diagrams of Fig. [4] were calculated by expanding the virtual momentum k\ in soft, 
hard and (anti-)collinear regions (the latter with respect to both legs p and p). This involves 
decomposing the virtual momentum in a Sudakov decomposition, according to 


K = \ ( n - ‘ fc i) n + + \ ( n + ‘ fc i) n - + K± > 


(4.24) 


where the vectors n± were chosen to be parallel to the p and p directions, respectively. In fact, 
they correspond to the vectors hi introduced here in Section [2] It was then found that the 
only contributions to the K-factor originated from the hard and (anti-)collinear regions. The 
connection with the present study is as follows. The contributions from the hard region consist 


of the first two terms on the right-hand side of Eq. (4.23), namely to the one-loop non-radiative 


amplitude dressed by an additional emission, and the derivative of the same amplitude. The 


remaining term in Eq. (4.23) comprises contributions from the collinear regions associated with 


each leg. We see this directly in Eq. ( 3.14| ), in which the terms on the right-hand side are explicitly 
proportional to an overall scale factor of (2 p-k)~ e or (2 p-k)~ e . This is the only scale that survives 
in the collinear region, and it is precisely these contributions that the original Low’s theorem [48] 
fails to capture, but which are discussed extensively in Ref. m (see also Ref. [ 561 ]). Note that 
the exact correspondence is achieved because we have chosen here the m vectors in accordance 
with a method-of-regions, or effective-theory approach. 

It is also interesting to examine our results in light of the recently proposed next-to-soft the¬ 
orems in gravity m and gauge theories |80| . discussed further in Refs. HEHMl MSI E2J. First, 
one may note that the next-to-soft theorems as presented in Ref. [791 [80] involve an explicit cou¬ 
pling to the total angular momentum of each outgoing leg, a structure which is not immediately 


apparent in Eq. (2.29). This is ultimately due, however, to the fact that Eq. (2.29) is presented 


in a form which is most straightforward for explicit calculation of the various contributions. Nev¬ 
ertheless, one does see this structure emerge at loop levej^] in Eq. (3.5), where the second and 
third terms on the right-hand side combine to give 




(4.25) 


(i) 

where is the spin angular momentum associated with leg i, defined in Eq. (|3.8[), and 


r <i) .1 d d 

L \,l = X wPiv ~ XivPi/j, = ~ Viv ~d^ 


(4.26) 


1 The tree-level case has already been examined in Ref. 
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is the orbital angular momentum. Eq. (4.25) is precisely the coupling to the total angular momen¬ 
tum observed in Ref. [80]. Note, however, that the final term in Eq. (3.5) involves a coupling to 


the spin angular momentum only, and thus corresponds to an explicit breaking of the next-to-soft 
theorem at loop level. As discussed above, this breaking is associated with collinear effects. 

A particular point of discussion in recent literature is whether or not the next-to-soft theo¬ 
rems of Refs. [791 80] receive corrections at loop level. This has been related to the sequential 
order in which one performs the dimensional regularization and soft expansions, with [85] ad¬ 
vocating performing the soft expansion first. This was further discussed in Refs. [82l 88l 189] , 
with Ref. [82] strongly arguing that the expansion in the dimensional regularization parameter e 
should be carried out first. A concrete example was examined in Ref. [63], namely the real-virtual 
corrections to the Drell-Yan K factor, that we have also considered here. There, it was found that 
contributions with a logarithmic dependence on the radiated gluon momentum k are generated 
in the amplitude from the collinear region. These are needed to reproduce the known K factor, 
and would be manifestly absent upon carrying out the soft expansion before the e expansion. In 
the present context, this corresponds to expanding in k before carrying out the integration over 
the virtual gluon momentum k \. The presence of such terms indicate that there are indeed loop 
corrections to next-to-soft theorems. 

The results of this paper, building upon the earlier work of Ref. [ 50] , provide for the first time 
a description of this breaking in terms of a universal factorized expression. It is then interesting 


to trace these loop corrections in the NLP expression for the amplitude, Eq. (2.29). As Ref. 


makes clear for the Drell-Yan example, the terms which are problematic from the point of view 
of the next-to-soft theorems, and which depend upon the ordering of the soft and e expansions, 
are those with an overall power of (p • k)~ £ , or (p ■ k)~ e . As already noted above, these stem 
from the collinear regions associated with the corresponding external legs, and appear here in 
the second line of Eq. (2.29), through the radiative jet function J^. This is not surprising: the 


failure of the next-to-soft theorems at loop level means that extra information must be inserted 
by hand in order to patch up the collinear regions. This explains the presence of the jet emission 
function, and its definition makes clear that such terms are universal. 


5 Conclusion 

In this paper, we have considered the generalization of the soft-collinear factorization theorem 
for gauge-theory scattering amplitudes to include corrections which are responsible for next-to- 
leading power threshold logarithms in high-energy cross sections. We concentrated on electroweak 
annihilation cross sections, using the Drell-Yan process as an example, since in this case final 
state radiation is forced to be (next-to-) soft near threshold, and the problem of NLP threshold 
logarithms can in principle be attacked using the tools provided by the Low-Burnett-Kroll-Del 
Duca theorem. The organization of these logarithms at finite order, and possibly to all orders in 
perturbation theory, is likely to play a crucial role in improving the precision of collider physics 
predictions, and it might also have a number of more formal applications, concerning the infrared 
limit of scattering amplitudes beyond divergent contributions. 

Building upon the work of Ref. [50] . we have constructed an explicit expression for the 
radiative scattering amplitude, in which a (next-to-) soft gluon is emitted into the final state of 


the Drell-Yan process. Our basic result is best summarized in Eq. (3.3), which expresses NLP 


corrections to the amplitude in terms of known universal quantities that appear in the soft- 
collinear factorization formula, as well as a new universal quantity, the radiative jet function , 
defined in terms of an auxiliary ‘factorization vector’ n in a manner similar to ordinary (non- 
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radiative) jets. This function first appeared in Ref. |50j . and is computed here for the first time 
at the one-loop level. We have also refined the treatment of Ref. [50] by introducing the explicit 
factorization of soft modes, including the correct treatment of the double counting of soft-collinear 
singularities, and by studying the role played by the factorization vectors to*. In particular, we 
have noted that it is possible, in this case, to work with light-like reference vectors, nf = 0, 
without introducing any further double countings of poles, or other inconsistencies, thanks to the 
simple renormalization properties of the proposed NLP expression for the amplitude, Eq. (3.3). 
This has led us to a considerably simplified formula, Eq. (3.5), in terms of light-like reference 
vectors, which can be directly applied to construct NLP-accurate cross sections. 

In order to test our generalized factorization formula, we have used it to reproduce the 
known NLP contributions in the abelian-like real-virtual contribution to the Drell-Yan K-factor 
at NNLO. Purely real emission (of two gluons) at this order was already understood for this 
process in Ref. 153, using results based only on the soft expansion, pushed to next-to-soft level. 
It is only when virtual gluons are present that one becomes sensitive to collinear effects. Thus, our 
exercise is a highly non-trivial check of Eq. (3.5), in which all parts of the formula are tested. We 
find that we do indeed reproduce all NLP threshold contributions, including even the constant 
term in the finite part, which is not logarithmically enhanced. This suggests that Eq. (2.29) can 
be used even at finite orders to obtain good approximate higher-order cross sections at NNLO and 
beyond. We emphasize that, while our current results were derived for the abelian color structure, 
using a QED-like expression for the current, the formalism can be naturally generalized to the 
full non-abelian theory, by using the non-abelian current in Eq. (2.16). When this is done, the 
process of Higgs production via gluon fusion will become accessible to our treatment, as well as 
other gluon-initiated cross sections. 

Our results are an essential first step to eventually develop a complete resummation formalism 
for NLP threshold logarithms, but, to this end, much work remains to be done. The next step 
will be to precisely work out and test the non-abelian version of the present formalism, which 
will immediately lead to an application to Higgs production in the gluon fusion channel. We will 
then need to consider in more detail the factorized expressions that arise at the level of cross 
sections, rather than amplitudes, including a fully systematic treatment of multi-gluon phase 
space. Finally, a non-trivial step will be the generalization to final state jets, which will require a 
slightly different treatment, due to the importance of hard-collinear gluon contributions, which 
cannot be reached through the soft expansion. Work on all these issues is in progress. 
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Figure 4: Feynman diagrams for the matrix element squared contributing to Drell-Yan production 
at NNLO, and involving one real and one virtual emission. Diagrams obtained by interchanging 
pop and/or complex conjugation are not shown. 
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A The abelian-like real-virtual NNLO Drell-Yan K-factor 


In Eq. (4.7), we quoted a result for the Cp part of the real-virtual contribution to the NNLO 


Drell-Yan A-factor. This has not previously appeared in the literature, so we believe it is useful 
to collect a few details regarding the calculation. In particular, we show here, in Fig. |4j the 
relevant Feynman diagrams for the squared matrix element, and we list below their individual 
contributions to the cross section at NLP in the threshold expansion. We have carried out the 
Dirac traces, the integration over the virtual momentum k\. and the Passarino-Veltman reduc¬ 
tion [MS using FeynCalc [108/ Standard results for one-loop scalar integrals may be taken, 


for example, from Ref. [ 109j . After substituting the phase space parametrization in Eq. (4.5), 


one may expand the squared amplitude in powers of (1 — z), where NLP accuracy corresponds 
to keeping the first subleading correction. One may then straightforwardly integrate over the 


real-gluon phase space, using Eq. (4.6). Discarding transcendental constants, as in the text, the 


results for individual Feynman diagrams at NLP accuracy are given by 
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(A.3) 

192, 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


where we have used the definitions of Eq. (4.8), and we have multiplied each diagram by four, 
in order to take into account complex conjugate diagrams, as well as those obtained from Fig. [f] 
by interchanging p and p. 

It is interesting to note that there are eikonal terms in graphs (e)-(h), and they cancel when 
all such contributions are added together. This had already been noted in Ref. |64j . and the fact 
that such terms appear in individual diagrams is an artifact of our (Feynman) gauge choice. 
Summing all diagrams together, one obtains the result of Eq. (4.7). Note also that we have not 
included UV counterterms in the above calculation: if included, these contributions cancel, in a 
manner reminiscent of the Ward identity that prevents QCD corrections from renormalising the 
(QED) photon-fermion vertex. 
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